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Abstract 

In the context of brane-world theories, the production cross-section for black holes may 
be greatly enhanced. Such black holes can in principle lead to detectable signals via 
their Hawking evaporation to brane-localized modes. We calculate, in the semiclassical 
approximation, the leading corrections to the energy spectrum (the greybody factors) 
for decay into scalar fields, as a function of the number of toroidally compactified extra 
dimensions, and partial wave number. 



1 Introduction 
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It is a remarkable fact that in brane- world theories the true scale, M*, of quantum gravity 
may be substantially lower than the traditional Planck scale, M p i, possibly approaching 
the TeV-scale, and this radical departure from the standard picture is not excluded H 
(for earlier works on brane theories see 0). This observation has naturally excited a 
large amount of activity investigating both the structure of these theories and their 
experimental signals j3| . One of the most striking consequences of lowering the Planck 
scale to the TeV region is that the properties of small black holes are substantially 
altered Q. A black hole of given mass M is now much lighter, larger, and colder than 
a usual black hole of the same mass, provided only that the Schwarschild radius rjj of 
the black hole is smaller than the size of the extra dimensions r# < R. In this limit, the 
black hole is well-described as a (4 + n)-dimensional black hole centered on the brane, 
but extending out into the n extra dimensions. The horizon radius r# of such a black 
hole is [§§, 

l/(n+l) 

M*\Mj \(n + 2)n( n+1 )/ 2 ) 

(Note that, following common practice, we work in the approximation that the brane 
tension itself does not strongly perturb the (4 + n)-dimensional black hole solutions.) 

In particular it is likely (though not proven) that black holes are much easier to 
produce, with production cross-section at parton-parton c.o.m. energies y/s close to the 
geometrical cross-section of a black hole of mass M = \fs || 0] (for supporting evidence 
see, e.g., ||, for claims to the contrary, see ||) 

, i (m /sr(^)l\ 2/( " +1> 

If this is the case, then for M* ~ 0{ TeV) there are striking consequences for the high- 
energy interactions of cosmic rays HlOU , and, moreover, the LHC will become a 'black- hole 
factory' |,0,[rj. 

After such black holes are produced they decay by the emission of Hawking radiation. 
It is expected that black holes produced by a collision on our brane, e.g. pp — > BH + X at 



a hadron collider, will decay mostly to particles on our brane |12| , and thus be indirectly 
observable via this characteristic Hawking radiation. This radiation is usually described 
as 'thermal' in character with a temperature 

Tbh = (3) 
Air th 

(in Gn = }vb = c = h = 1 units). However, because of the non-trivial metric in the region 
exterior to the horizon there exists an effective potential barrier in this exterior region. 
This potential barrier backscatters a part of the outgoing radiation back into the black 
hole, the amount depending on the energy of the radiation. Thus the original blackbody 
radiation is modified by a frequency-dependent filtering function, <r{uS), caused by the 
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gravitational potential of the black hole. The function, <r(ou), is known as the 'greybody 
factor'. The black hole differential decay rate into particles of energy uj is then given by 
the Hawking formula |13| 



dE{u) ^ , , uj d n+3 k 

£,6 



where i labels the angular momentum quantum number, b labels any other quantum 
numbers of the emitted particle, as well as the particle type, and in the phase-space 
integral \k\ — u for a massless particle.[| 

Greybody factors are important theoretically because they encode information on 
the near horizon structure of black holes. Indeed one of the most exciting features of 
BH production at the LHC would be the opportunity to investigate Hawking emission 
in detail for clues as to how the infamous apparent violation of the laws of quantum 
mechanics by black hole evaporation (the 'information paradox' of black holes ||14j| ) is 
resolved - if it is, that is! 

Greybody factors can be important experimentally because they modify the spectrum 
in the region where most particles are produced thus altering the characteristic spectrum 
by which we hope to identify a 'BH event'. In particular the functional dependence of 
<Jb{^) on the energy uj depends on the spin of the emitted particle, and on whether it is 
brane-localised or free to propagate in the bulk of the extra dimensions. 

One can compute the greybody factor by first computing the absorption cross section 
for the appropriate type of particle incident on the background metric that describes the 
brane black hole. This is because the greybody factor in the Hawking formula for the 
emission rate of a given type of outgoing particle, b, at energy uj equals the absorption 
cross section <Tb,abs(w) for a particle of type b incoming at energy uj. In fact it is this 
property which implies the greybody factors do not invalidate the thermal nature of 
the black hole. Since the outgoing transmission and incoming absorption coefficients are 
equal to one another, equilibrium still occurs if the black hole is placed in a heat bath. 

The semiclassical calculation of Hawking emission is only reliable when the energy of 
the emitted particle is small compared to the black hole mass uj <C M, since only in this 
case is it correct to neglect the back reaction of the metric during the emission process. 
This in turn requires that the Hawking temperature Tbh ^ M, which is equivalent to 
demanding that the black hole mass M 3> M*, as can be seen from Eqs.([|) and (|l]). 
Inevitably this condition breaks down during the final stages of the decay process, but 
for those black holes of initial mass larger than M* most of the evaporation process is 
well- described by the semi-classical calculation. 

In this paper we calculate the greybody factor for the simplest case of scalar particles, 
both free to propagate in the bulk, and localized on the brane (e.g., the Higgs). Section 2 
contains a discussion of the metric used in calculating the greybody factors and the 
quality of the approximations used. To orient the reader through our analysis, Section 3 



*An alternate form for the decay rate that is sometimes useful in the massless particle case in- 
volves the absorption probability |^(tj)| 2 whose relation to <j(u>) is given in Eq. (|l3|) : dE(uj)/dt = 
m b \ A(lu)\ 2 ujduj / [2ir(exp (oj /Tbh) T !)]• Here m is the azimuthal quantum number. 
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presents the calculation of bulk scalar emission in (4 + n)-dimensions, in the simple case 
of S-wave emission. Section 4 performs the calculation of the greybody factor for bulk 
scalar emission for arbitrary partial wave. The primary result of this section, the general 
formula for the greybody factor, is given in Eq.(^). Section 5 contains the calculation 
of brane-localized scalar emission, where the brane is embedded in a (4 + n) -dimensional 
bulk. The primary results of this section are given in Eqs.([|7]) and d48f) . We conclude in 
Section 6. For ease of use, the results of our calculations for the greybody factors, for 
the most important (I = 0, 1, 2) angular momentum modes, and for the n = 2, 4, 6 extra 
dimensions are collected in Tables 1 and 2. 



2 Black hole metrics and greybody factors 

Let V = (2nR) n be the volume of the n extra dimensions, which are here taken to 
be of common radius R. Then Gauss' Law relates the 4d Planck mass M p i to the 
new fundamental scale of gravity by M 2 X = V n M 2+n . A black hole of mass M <C 
M pl (M pl /M Sf ) ( - 2+n ^ n ~ 10 15+32/n GeV is, for distances r < R, well approximated || 
by the (4 + n)-dimensional Schwarzschild black hole with line-element 

ds 2 = -h(r)dt 2 + h(r)~ l dr 2 + r 2 dQ 2 2+n , (5) 

where 

h(r) = I ~ (-Y + \ (6) 



r 

with the horizon radius, r#, given in Eq.(|l|). In Eq.(^) the angular part is 

dQ 2 2+n = d9 2 n+1 + sin 2 9 n+1 {de 2 n + sm 2 9 n ( ... + sin 2 9 2 {d9{ + sin 2 6 1 d<p 2 ) ... ) V (7) 

with < (p < 2tt and < Q{ < 7T, for i — 1, n + 1. However, because of the compacti- 
fication of the extra dimensions, the metric for this black hole at distances r 3> R goes 
over to that of the usual 4-dimensional Schwarzschild solution 

-l 



d? = -[l-™)dt*+(l-^-) dr 2 + r 2 dn 2 2 



The matching of the two expressions, Eqs.(^) and (Q) at r ~ R is of course just 
an approximation to the exact metric for such a black hole. An exact expression is not 
necessary to derive the form of the greybody factor to a very good approximation for 
the energies of interest. 

To understand this consider the case where a scalar field propagates in the full (4+n)- 
dimensional Schwarzschild black hole background. The scalar wave equation g IJ (pj-j = 
in this background is separable if we make the ansatz 

<f>(t,r,ei,(p) = e^R^(r)Y e (n), (9) 

where i^(fi) is the (3 + n)-spatial-dimensional generalisation of the usual spherical 
harmonic functions depending on the angular coordinates [15]. Upon substituting this 
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ansatz, the scalar wave equation implies a second-order differential equation for the radial 
wavefunction R^Ar): 



h(r) d 

r n+2 ^ r 



dr 



+ 



u 2 -^£(£ + n + l) 



R = 0. (10) 



The greybody factor for BH decay into scalars will be found from the solutions to this 
equation. 

We can transform this equation to a more convenient form by defining a new ( "tor- 
toise" ) radial coordinate y by 

ln/i(r) dy 1 

r H +1 (n + 1 ) dr ~ K r ) rT 
in terms of which the radial equation becomes 

2 £(£ + n + l)h(r) 



- v ; =± zE. — (]]) 
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R(y) = 0. (12) 



Since the coordinate position r# of the horizon is defined by the (largest) solution to 
h(rn) = 0, the horizon in terms of y is at y — > — oo, as can been seen from Eq . (p|) . 
Eq.(|12"D is analogous to the Schrodinger equation for a particle in an effective potential. 



Non-trivial backscattering occurs in this metric when the 2nd of the two terms in the 
[] -parentheses in Eq. (|12"|) is comparable or larger than the 1st. For black holes at the LHC 
the typical energies of emitted radiation we are interested in range from u> max ~ Tbh = 
(1 + n)/47rrff ~ (few) 100 GeV to a minimum of ~ (few) GeV (this minimum might be 
set by the energy threshold of the detectors - the precise value will not matter) . For the 
range of parameters, n, and M*, of interest, this gives a range of u from u> max < l/2r H to 
Umin ~ l/200r#. Thus, inspecting the potential terms of Eq.fll^) we see that the potential 
is only large enough to lead to backscattering over the range of distances r from (9(l)r# 
to O(100)r H . 

A similar analysis can be performed for the 4d asymptotic metric of Eq.(||). In this 
case one finds that significant backscattering of the quanta of interest would only occur 
when distances were of order r ~ M/M^ C C fl. But at such distances the 4d 
asymptotic form of the metric is not applicable, being replaced by the (4 +n) -dimensional 
metric used above. Instead what happens is that the change-over from the (4 + n)- 
dimensional to 4-dimensional regime, at distances of order r ~ R, only significantly 
backscatters quanta of energy 10~ 3 eV or less. In other words, the backscattering due to 
the effect of the compactification of the extra dimensions only affects very low energy 
quanta which are experimentally irrelevant. Even more so this applies to the change in 
backscattering at distances r ~ R due to the difference between the approximate form 
of the black hole metric, Eqs. ([|) and (|^), and the exact black hole solution. 

Overall, the lesson is that upto corrections that only apply for very low energy, 
the greybody factors for brane black holes may be calculated using purely the (4 + n)- 
dimensional regime of the black hole metric. In the remainder of this paper we will follow 
this procedure. Thus our task is now to return to the radial equation Eq. ( |j"U| ) or ( |T2"D and 
compute the greybody factor. 



4 



In particular, via the equivalence of <j(lj) to the absorption cross section, the greybody 
factor can be computed by first evaluating the absorption probability, |^4(a;)| 2 , from the 
ratio of the in-going flux at the future horizon to the incoming flux from past infinity 
(with the boundary condition that there is no outgoing flux at the horizon), and then 
using the generalised (4 + n) -dimensional optical theorem relation |TjJ 

2^+i)/2r> + l)/2] (2l + n + l)(l + n)! , „ 2 

= J\ W ( 13 ) 

between the absorption cross section a^(u) and the absorption probability \A\ 2 for the 
£'th partial wave. This formula includes a summation over the multiplicity of individual 
'azimuthal' components for each given partial wave £ (i.e. 2£ + 1 in 3 spatial dimensions, 
(£ + l) 2 in 4 spatial dimensions, etc.). A simple summation is appropriate because we 
deal with a non-rotating black hole. 

The radial equation, however, is not in general exactly soluble, therefore, we will em- 
ploy an approximation method based on splitting the radial domain into 'near- horizon' 
(NH) and 'far-field' (FF) regions. The solutions (satisfying appropriate boundary condi- 
tions) in these two regions will be computed and then matched in a transition region to 
find the complete solution. This procedure leads to an expression for |^(a;)| 2 correct in 
leading order in an expansion in the dimensionless quantity ujth- 



3 Bulk scalar emission: S-wave example 

Because the radius and temperature of the (4 + n)-dimensional black hole are always 
comparable, Eq.(|), the dominant scalar decay mode is that into the S-wave, £ = 0. We 
will first solve the problem in this case, which will also serve as an illuminating example 
for the full case studied in Sections 4 and 5. 

We will first compute the solution in the 'near-horizon' (NH) region. The radial 
equation for the S-wave, in terms of the y-coordinate defined in Eq. (Jl~T|), becomes 

(^- 2 +" 2 r 2(n+2) )R(y) = o. (14) 

By expanding close to the horizon r = rjj + 5r (5r <C r#), we obtain 5r ~ 
rf/exp f(n + l)r^ +1 ?/) as y — > — oo. Thus the radial equation Eq. flliD becomes in the 
near-horizon limit 

(|,+-vr 2) )^)=o, as) 

up to exponentially small corrections in y. The general near-horizon solution is therefore 

R NH (y) = A + exp(irl +2 Luy) + A_ exp(-ir™ +2 uy). (16) 

In order to calculate the greybody factor, we must impose the boundary condition that 
near the horizon the solution is pure in-going. Therefore, we need to set A + = 0. 
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We now turn to the far-field region which is defined by r ^> r#. In this limit, h(r) ~ 1 
and, by setting R(r) = /( r )/ r («+i)/2 ; Eq.(0) can be rewritten as 



d l f 1 df 
— - H + 

dr 2 r dr 



u 2 - 



(71 + 1) 



/ = 0, (17) 



which has the form of a Bessel differential equation WJ§. The general solution for the 
radial function i?(r) is therefore given by 

B B 

R FF{r) = J (n+1)/2 (u;r) + F ( „ +1)/2 (a;r) , (18) 

where Jf n+ iy 2 (ujr) and Y( n+ iy 2 (u;r) are the Bessel functions of the first and second kind, 
respectively 

As we will soon see, the two coefficients B + and B_, and more specifically their ratio, 
will help us define the greybody factor. To compute this ratio, we need to match the 
far- field solution, Eq. (|T8"D , on to the near- horizon solution, Eq.flTEp, in the intermediate 
region. To this end, we expand the near-horizon solution, in the regime ur <C 1 and 
r 3> r H, leading to 

( ur {TH\ n+2 \ [ . ur frH^ n+2 ^ 



^"M^-^twlfJ )^ A -\ l+ '—Af) }■ < 19 > 

We also expand the far-field solution, Eq. (|T8"D , in the same regime ur <C 1, which gives 



Bm /u\ {n+l ^ 2 B /2\ { - n+l)/2 T{^ ±1 ) 

"»M°Tm{i) (20) 

Matching the above expression with Eq. ([T9|) , we find the result 

¥± = j - { ~ ] (21) 

The reflection coefficient 1Z for scattering in the gravitational potential of the black 
hole Eq.(|5|) is defined as the ratio of the outgoing and incoming amplitude at infinity. 
To compute this, we expand the far-field solution Eq. flTjD , in the limit ur — » oo, which 
yields 



{B + - iBJ) (B + + iB_) — il lot— (n.+i)7r/4 

J~iFF{ r j — — / e V > H , e V / 

V2Tnor n + 2 V2iru;r n + 2 

and which, in turn, leads to the following expression for the reflection coefficient 



(22) 



outgoing amplitude B + — iB_ 
incoming amplitude B + + %B_ ' 



2 This expansion allows u < 1/rsr, the typical emitted energy, when M 3> M*, as is anyway required 
for the reliability of the semiclassical approach. 
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up to a purely imaginary phase that will drop out when the magnitude of 1Z will be 
computed. The absorption probability is then defined as 



W 2 -(1-| K |^|^, (24) 

where, in the final expression, we have expanded to leading order in (urn) • This is the fi- 
nal result for the S-wave greybody factor for scalars in a (4+n)-dimensional Schwarzschild 
black hole background. We may easily check that, for n = 1, we correctly reproduce the 
result 

\Af = \u 2 r H = ^-A H , (25) 



where Ah = 2ir 2 r H is the area of the horizon, presented in Ref. WE 



4 Bulk Scalar emission for £ > 

We will now generalise the above analysis in the case where the scalar modes are not 
spherically symmetric, i 7^ 0. For readers more interested in our final results, rather than 
the techniques used to solve the problem, we suggest jumping to Eq.(|57|) and following 
discussion. 

We will start with the derivation of the solution in the near-horizon zone. Starting 
from Eq. (|10D , and making a change of variable, we may write the scalar field equation 



in the form (here, we adopt the method of Ref. [[19 



, s d 2 R , r ,dR 



(lot) 2 i (£ + 1 + n) 



(n + l) 2 h(l-h) (n + l) 2 (l-h) 



R = 0. (26) 



Near the horizon, r ~ r# and the quantity (ur) 2 can be set equal to (uth) 2 - Then, by 
using the redefinition R(h) = h a (l — /i) /3 F(/z), the above equation takes the form of the 
hypergeometric equation 

d 2 F dF 
h(l - h) — + [c - (1 + a + b) h] — - abF = , (27) 
dli z dn 

with a = b = a + j3 and c = 1 + 2a, where 



— Ifr- * = 5 ^c+^r -<-»>'■ (28) 

Equation (p7|) has as a solution the hypergeometric function F(a,b,c; h). The criterion 
for the convergence of the hypergeometric function demands that TZ(c — a — b) > 0, which 



forces us to choose (3 — Then, the general solution of Eq . (|2q) may be written as [17] 



R NH (h) = A_h a± (1-hf F(a,b,c; h) + A + h~ a± (1 -hf F(a- c+ 1, b- c+ 1, 2 - c; h) . 

(29) 
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Expanding the above solution in the limit r or /i —> 0, and choosing a = a_, we 

obtain the result 



R 



NH 



(30) 



which again imposes the condition^ A + = 0. 

The derivation of the far-field-zone solution closely follows the analysis performed in 
the case where i = 0. The same redefinition of the radial function leads again to a Bessel 
equation whose general solution is now given by 



B B 



(31) 



Both solutions, near-horizon and far-field, need to be "stretched" and matched in 
the intermediate region. Expanding first Eq . (|3l|) , in the limit ur <C 1, gives 



R FF {r) 



B: r l 



w \fr(»+l)/2 



t+n+l 



r 



2 y+(n+l)/2 r ^ + n±l) 



UJ 



7T 



(32) 



The near-horizon solution needs to be "shifted" first and expressed in terms of 1 — h, 



before being expanded in the limit r ^> r#. By using a standard formula [17], we write 



R NH {h) = A^h c 



+ (i-/0 



1 T(l + 2a) r(l - 2/3) 



r(i + «-/3) 2 



F(a,b,a + b - c + 1; 1 - /i) 



1 _ /J r(l + 2a)r(2/3-l) 



F(c — a, c — 6, c — a — b + l;l — /i) 



1> + /?) 2 

We can now expand the above expression in the limit ft — > 1 and take 



.(33) 



p ( h \„A r y T(l + 2a) r(l - 2/9) , , ^r^+n+i T(l + 2a) r(2/? - 1) 

« -^-( — ) — 7^77- ^2 + A ~{ — ) 77 — — ^ • 34 



Matching the two solutions Eq.(|3"2"D and Eq. fl3~4"|) , we obtain the ratio 

2 x 2£+n+i T(£ + ^±±) 2 (e + 2±1) T(l - 2/3) r(a + /3) 2 



7rr(l + a-/3) 2 r(2/3- 1) 



(35) 



The definition of the reflection coefficient 1Z is still given by Eq . (|23|) . In turn, the ab- 
sorption probability can be written, in terms of B = B + /B_, as 



:i - \n\ 



2i (B* - B) 



BB* + i(B* - B) + 1 



(36) 



3 Note that the choice a = a + would have led again to Eq.(p0|) with A_ <-» A + , and therefore to the 
choice A_ = 0. As both values of a appear in the general solution Eq.(p9|), it is only a matter of choice 
which one of the two terms will be associated with the incoming mode. 
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Due to the fact that the argument of the Gamma functions appearing in the expression 
of 1Z are non-trivial complex numbers, we can not write the absorption coefficient in a 
simple way in the general case. However, the above expression can be further simplified 
in the limit loth <C 1, in which case BB* ^> i(B* — B) ^> 1, and, therefore, \A\ 2 may be 
written as 

4vr 2 /ujru\ 2e + n + 2 T(l + — ) 2 

\A\ 2 = f Wg l 1 I w (17) 

11 2*+dU j r(i + 4 T ) 2 r(£ + ^) 2 ' 1 J 

This is our major result for the case of bulk scalar fields. Eq.(^) nicely displays the 
leading functional dependence of the greybody factor on ujth for varying partial wave, £, 
and number of extra dimensions, n. For the case of £ = 0, this may be further evaluated 
to give 

2 fur H \ n+2 47T 
l^=o = {—) r[(n + 3)/2 ]2 > ^ 



in complete agreement with the result, Eq.(E¥), of our earlier analysis in the case £ = 0. 
In fact, for an s-wave massless bulk scalar, the absorption probability \A\ 2 has the exact 
form that allows the greybody factor cr(u) to reduce to the horizon area of the black hole 
in agreement with previous work p0j. The numerical value of the results for the cases 



£ = 0, 1, 2 for n = 2, 4 and 6 are shown in Table 1. 

If we fix the number of extra dimensions and vary only the angular momentum 
number, the absorption probability decreases as £ increases. This decrease is caused by the 
fact that the power of the expansion parameter (ur H) 2l+n+2 , in the leading term, increases 
with £. Since (our H ) <C 1, this means that \A\ 2 becomes more and more suppressed as 
£ increases. The numerical coefficient in front of the leading term also decreases with £ 
(see Table 1). The same behaviour is observed if we fix instead £ and vary n. 



Table 1 : 


Absorption probabilities for a (4 + n) bulk scalar field 


n = 2 


£ = 


\A\ 2 ~ l(u;r H y + . . . 




£= 1 
£ = 2 


| ,12 _ 2^/3 r(4/3)» / )6 

\^\ - (15) 2 r(5/6) 2 n \ Ur H) "I"--- 

| ,12 _ 2-2/3 r(5/3) 2 / )S , 
l^l — (105) 2 T(7/6) 2 n \ U 'H) -T ■■■ 


n = 4 


£ = 


\A\ 2 ^^(ur H f + ... 




£= 1 
£ = 2 


I ,12 _ 2«/« r(6/5) 2 / )8 , 

l^l — (105) 2 r(7/io) 2 n \ ur H) -I- • • • 

1 4 12 _ 2 2 /5 r(7/5) 2 / uo , 
l^l — (105) 2 3 4 r(9/10) 2 71 \ U 'H) I"--- 


n = 6 


£ = 


I^l 2 ^ (10 4 5)2 M 8 + ... 




£= 1 
£ = 2 


I ,12 223A r(8/7) 2 / UO . 
I^l — (105) 2 3 4 T(9/14) 2 n K^'H) -I"--- 

| ,12 _ 2^ r(9/7) 2 / U 2 , 

l^l — (ii55) 2 3 4 r(n/i4) 2 71 y u 'H) -r ■■■ 



9 



5 Brane-localized scalar emission for I > 



We now turn to the study of the case where the scalar field is confined on a four- 
dimensional brane embedded in a (4 + <i)-dimensional Schwarzschild spacetime. The 
scalar field propagates in a four-dimensional background whose metric tensor is given 
by the induced metric at the location of the brane. The induced metric follows from the 
(4 + (i)-dimensional one by fixing the values of the extra angular coordinates: 8 n = tt/2 
for n > 2, and it may be written as 

ds 2 = -h{r)dt 2 + h(r)- l dr 2 + r 2 {d8 2 + sin 2 9 dtp 2 ), (39) 

where h(r) is still given by Eq.(|6]). The scalar field equation may be separated in the 
same way 

<P(t,r,9,<p) = e- ta *R ut (r)Y e (n), (40) 

where Ye(Q) are now the usual 3-dimensional spherical harmonic functions. The above 
ansatz allows us to write the equation for the radial part as 



h(r) 



d 
dr 



2 dR 
hyr) r — — 
dr 



+ 



h{r) 



£(£+!) 



R = Q. 



(41) 



The presence of the metric function h(r) makes once again the derivation of the general 
solution extremely difficult. We will follow the same method as in the previous section 
and compute the solution in the two radial domains, near-horizon and far-field, which 
will then be "stretched" and matched in the intermediate region. 

Having become familiar with the analysis, we will proceed to derive directly the 
solution in the general case £ > 0. We start with the solution in the 'near- horizon' (NH) 
region. In terms of h, the radial differential equation now takes the form 



h(l- 



1N d 2 R 



(2n + l) 
(n + l) 



h 



dR 
dh 



+ 



(ur H f 



£(£+1) 



(n + l) 2 h(l-h) (n + l) 2 (l-/i) 



R = 0. (42) 



By using the same redefinition R(h) = h a (l — /i) /3 F(/i), the above equation assumes the 



standard form of a hypergeometric equation with indices a 
and c = 1 + 2a, where 



a + /3 + 



(n+l) ' 



a± = ± 



n + l 



01 



1 ± J (21 + l) 2 - A(ur H )' 



a + (3 



(43) 



2(n + l) 

The criterion for the convergence of the hypergeometric function, lZ(c — a — b) > 0, 
demands again f3 = The general solution of Eq.(|42]) has again the form of Eq. (f29|) . 
Expanding near the horizon and imposing the condition that only incoming waves exist 
near r ~ r#, we find that A + = for a = a_. At this point, we can also "shift" the 
solution and write it in terms of 1 — h, in the same way as in the previous section. If we 
finally expand for r ^> r#, or equivalently h — > 1, we obtain the solution 

r(l + 2a) r(l - 2(3 - ^j) 



R NH (h)~A„ — 



/r^y+i r(l + 2a)r(2/? + ^-l) 
V r ) T(a + (3+^ l )Y{a + (3) 



(44) 
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The far-field-zone solution can be easily found to be given in terms of the Bessel 
functions J^ + i/ 2 (co>r) and Fg +1 / 2 (W). Expanding the general solution in the limit wr < 1, 
finally gives 



Rff(t) 



i+l/2 



■H'+l 



2^+1/2 r(^+|; 



71 



(45) 



Matching the two asymptotic solutions, we obtain the ratio 



tor H 



71 



r(l + a-/3)r(l + a-/3-^)r(2/3+^-l 



which can be used to determine the absorption coefficient according to Eq. 
however, obtain a simplified expression, in the limit uoru <C 1, which reads 



, (46) 
We may, 

(47) 



(n + i) 2 V 2 ; r(i + £)2r(i + f±i)2- 

The expression of the absorption coefficient \A\ 2 for the values n = 2,4 and 6 of the 
number of extra dimensions, and £ = 0, 1, 2 of the angular momentum number, are 
shown in Table 2. 

Finally, employing the relation between the absorption probability and the greybody 
factor, Eq. (p~3|) , leads to 

4^(21 + 1) r(^) 2 r(l + ^ (ur B \», 

r H • 



(n + iy r(i + ^r(i + f±i)2 V 2 J 



(48) 



In this equation, we have set n = in Eq.(|T^) as the partial waves are purely confined 
to the 3-spatial-dimensional brane, the only dependence on n being in \A\ 2 . 



Table 2 : 


Absorption probabilities for a (4D) brane scalar field 


n = 2 


£ = 




~ 4 (ujr H ) 2 + ■ ■ ■ 




£ = l 




^ff M 4 + ... 




i = 2 




^ (1 5 )2 (Wff) 6 + . . . 


n = 4 


£ = 




~ 4 (oor H ) 2 + . . . 




£ = l 




_ 4 r(2/5) 2 r(6/5) 2 , )4 , 

- 25 T(8/5) 2 i^H) + ■■■ 




1 = 2 




* + . . . 


n = 6 


£ = 




~ 4 (uor H ) 2 + ... 




£ = \ 




~. 4 r(2/7)2r(8/7)2 , u , 

- 49 r(io/7) 2 \ ur H) -r ■■■ 




£ = 2 




~ 4 r(3/7) 2 r(9/7) 2 / )6 

— (21) 2 r(12/7) 2 \ u ' H) 1- • • • 
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Keeping n fixed and varying i, we see that the absorption probability decreases once 
again as t increases (see Table 2): the dominant term becomes more and more suppressed 
by extra powers of (urn) and its numerical coefficient also decreases (the introduction of 
the multiplicity of states with the same angular momentum number i does not change 
this behaviour). If we fix instead t and vary n, a radically different behaviour, from 
the one observed in the case with a bulk scalar field, emerges: the leading term, in 
the expansion of \A\ 2 in powers of (urjj), remains the same, since it is n- independent, 
while its coefficient increases as n increases. In other words, for a given partial wave, 
the absorption probability, and therefore the greybody factor, increases as the number 
of extra dimensions being projected on the 3-brane also increases. 

6 Conclusions 

In this paper, we have studied the problem of scalar emission in a spherically symmetric 
.D-dimensional Schwarzschild black hole background. The cases of the emission of a 
(4 + n)-dimensional bulk scalar field and of a 4-dimensional brane-localized scalar field 
were studied separately and the greybody factor was determined in each case. This 
quantity causes the spectrum of Hawking radiation to deviate from the black body 
spectrum as it strongly depends on the energy of the particle mode emitted. Moreover, it 
encodes information about the gravitational background and thus on the number of extra 
dimensions that might exist (both in the case where the scalar field is free to propagate 
in the (4 + n)-dimensional bulk or when the field 'feels' the existence of extra dimension 
only through the induced metric on our 3-brane). 

The differential equation for the radial part of the scalar field can be solved by an 
approximation method valid in leading order in (curjj): the solution of this equation was 
found in the 'near-horizon' and 'far-field' region and were subsequently 'stretched' and 
matched in an intermediate regime. This matching allows us to determine the absorption 
coefficient for scattering in the black-hole background, which then leads to the greybody 
factor, cr(cj), according to Eq. (p~3|). 

We first focused on the case of a bulk scalar field propagating in a (4+n)-dimensional 
Schwarzschild black hole background. The general form for the amplitude probability was 
determined and an analytical, simplified version was also presented that allowed us to 
display the leading functional dependence on the expansion parameter (ooth) in terms 
of the number n of extra dimension and the angular momentum number £. Our results 
in this case are presented in Eq.(^) and in Table 1. 

The most phenomenologically interesting case is the one of a scalar field that is con- 
fined on a 3-brane and propagates in the induced spacetime background of a black hole 
(which is necessarily higher-dimensional). The functional form of the resulting greybody 
factor depends only on the angular momentum number, £, through (co>r#) 2f+2 . The de- 
pendence on the number of extra dimensions is entirely contained within the coefficient 
of this leading term. Our primary results are given in Eq.(f48|) and in Table 2. 

It is tempting to compare the results derived in the case of a brane-localized scalar 
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field (which nevertheless is part of a higher-dimensional manifold) with those valid in 
the case where a purely 4-dimensional scalar field propagates in a Schwarzschild black- 
hole background. Both cases lead to the same value of \A\ 2 for an s-wave and, thus, no 
distinction can be made between the two backgrounds. However, for higher partial waves, 
the value of the absorption probability in the case of a brane scalar field is always larger 
than the one for a purely 4D field, a fact which in principle can be used to distinguish 
between the two cases. 



In a companion paper we employ the techniques developed in the current paper 
to derive the greybody factors for higher-spin fields localized on a brane. This allows 
allows us to discuss the physics and phenomenology of black hole decay, as might be 
observed at the LHC. 



Appendix 

For completeness we here present some of the calculation of scalar emission in the case 
where the scalar field propagates in a purely 4-dimensional Schwarzschild background 
without the assumption of the presence of extra dimensions. The expression of the absorp- 
tion coefficient \A\ 2 can be easily found by first putting n = (the number of projected 
extra dimensions on the 4-dimensional plane) in the result for the ratio B + /B_. Then, 
Eq.(BBF) becomes 



B_ 



/ 2 \2wr(!+|f(f + |)r(i-2/?)r(a + ^ 



\ur H 



7rr(l + a-/3) 2 r(2/3-l) 



where now 



a = —Horn , 



l-J(2l + iy-A(ujr H y 



(49) 



(50) 



By using then Eq.(^), we obtain the result for \A\ 2 which, in its simplified form, reads 



\Af 



21+2 



r(£+i) 



r[i + 



r[2£ + 2] 2 



(51) 



Note that the same results follow by putting n = in all the expressions of section 2.1, 
as expected. 

We display the results for the absorption coefficient in the pure 4-dimensional case 
and for the values £ = 0, 1, 2 of the angular momentum number in Table 3. 



Table 3 : Absorption probabilities for a (4D) scalar field 


£ = 
£=1 
£ = 2 


\A\ 2 ~4(tur H ) 2 + ... 
\A\ 2 ~ | (ujr H y + ... 
\A\ 2 ^-L, K) 6 + . . . 
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